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The production of final state photons in hadronic Z-boson decays can be used to 
i-£h \ study the quark-to-photon fragmentation function D q ^ 1 (z, /xf). Currently, two dif- 

i ' ferent observables are used at LEP to probe this function: the 'photon' + 1 jet rate 

and the inclusive photon energy distribution. We outline the results of a calculation of 
the 'photon' + 1 jet rate at fixed 0(aa s ), which yield a next-to-leading order deter- 
mination of the quark-to-photon fragmentation function D q ^-y(z, ^p). The resulting 
predictions for the isolated photon rate and the inclusive photon spectrum at the same, 
fixed order, are found to be in good agreement with experimental data. Furthermore, 
we outline the main features of conventional approaches using parameterizations of 
the resummed solutions of the evolution equation and point out deficiencies of these 
currently available parameterizations in the large z-region. We finally demonstrate 
that the ALEPH data on the 'photon' + 1 jet rate are able to discriminate between 
different parameterizations of the quark-to-photon fragmentation function, which are 
equally allowed by the OPAL photon energy distribution data. 
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1 Introduction 



The production of final state photons at large transverse momenta is one of the key ob- 
servables studied in hadronic collisions. Data on high-p-r photon production yield valuable 
information on the gluon distribution in the proton, while the presence of photons in the 
final state represents an important background source in many searches for new physics. A 
good understanding of direct photon production within the context of the Standard model 
is therefore essential. 

Photons produced in hadronic collisions arise essentially from two different sources: the 
direct production of a photon off a primary parton or through the fragmentation of a hadronic 
jet into a single photon carrying a large fraction of the jet energy. The former gives rise to 
perturbatively calculable short-distance contributions whereas the latter is primarily a long 
distance process which cannot be calculated within perturbative QCD. It is described by 
the process-independent parton-to-photon fragmentation function |L| which must be deter- 
mined from experimental data. Its evolution with the factorization scale fip can however be 
determined by perturbative methods. 

Directly emitted photons are usually well separated from all hadron jets produced in a 
particular event, while photons originating from fragmentation processes are primarily to be 
found inside hadronic jets. Consequently, by imposing some isolation criterion on the photon, 
one is in principle able to suppress (but not to eliminate) the fragmentation contribution to 
final state photon cross sections, and thus to define isolated photons. However, recent analy- 
ses of the production of isolated photons in electron-positron and proton-antiproton collisions 
have shown that the application of a geometrical isolation cone surrounding the photon does 
not lead to a reasonable agreement between theoretical prediction and experimental data. 

An alternative approach to study final state photons produced in a hadronic environment 
is obtained by applying the so-called democratic clustering procedure 0. In this approach, 
the photon is treated like any other hadron and is clustered simultaneously with the other 
hadrons into jets. Consequently, one of the jets in the final state contains a photon and is 
labelled 'photon jet' if the fraction of electromagnetic energy within the jet is sufficiently 
large, 

z - -p p > Zcat, UJ 

-£->EM + &HAD 

with z cnt determined by the experimental conditions. This photon is called isolated if it car- 
ries more than a certain fraction, typically 95%, of the jet energy and said to be non-isolated 
otherwise. Note that this separation is made by studying the experimental z distribution 
and is usually such that hadronisation effects, which tend to reduce z, are minimized. 

This democratic procedure has been applied by the ALEPH collaboration at CERN in an 
analysis of two jet events in electron-positron annihilation in which one of the jets contains a 
highly energetic photon [[| . In this analysis, ALEPH made a leading order determination of 
the quark-to-photon fragmentation function by comparing the photon + 1 jet rate calculated 
up to O(o) M with the data. The theoretical basis on which the measurement of the 'photon' 
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+ 1 jet rate relies, is an explicit counting of powers of the strong coupling a s in both the 
direct and the fragmentation contributions, no resummation of In fj/p is performed. We shall 
refer to this theoretical framework as the fixed order approach. In Section 0, we will describe 
the main features of the leading and next-to-leading order calculation of the photon + 1 jet 
rate following this fixed order approach and see how the obtained predictions compare with 
the available data. 

More recently, the OPAL collaboration has measured the inclusive photon distribution 
for final state photons with energies as small as 10 GeV Q. This corresponds to the photon 
carrying a fraction of the beam momentum, x 7 , to be as low as 0.2. They have compared 
their results with the two model-dependent predictions of GRV []|] and BFG f| and found 
a reasonable agreement in both cases when choosing the factorization scale [ip = Mz- In 
Section |2| we shall compare the predictions obtained for the inclusive rate within our fixed 
order approach with the OPAL data too. 

The model predictions [||, |6| mentioned above are based on a resummation of the loga- 
rithms of the factorization scale fip and naturally associate an inverse power of a s with all 
fragmentation contributions. We shall refer to this resummation procedure as the conven- 
tional approach. In Section |3| we shall present the main features of this approach, outline the 
results obtained for the photon + 1 jet rate in this approach while using either the GRV or 
BFG parameterizations for the photon fragmentation function and show how these compare 
with the ALEPH data. Finally Section |4] contains our conclusions. 



2 The photon + 1 jet rate in the fixed order framework 
2.1 The photon + 1 jet rate at 0{a) 

In the fixed order framework, the cross section for the production of isolated photons receives 
sizeable contributions from both direct photon and fragmentation processes. More precisely, 
the distribution of electromagnetic energy within the photon jet of photon + 1 jet events, 
for a single quark of charge e q , at 0(a) in the MS-scheme, can be written as ||, 

+ R A 5(1 - z) + ..., 

(2) 

is the coefficient function 
corresponding to the lowest order process e + e~ — > qqj. It is defined after the leading 
quark-photon singularity has been subtracted and factorized in the bare quark-to-photon 
fragmentation function in the MS scheme. The non-perturbative fragmentation function is 
an exact solution at 0(a) of the evolution equation in the factorization scale /j,p, 

D q ^(z,fz F ) = °^P^{z) log (41 + D q ^(z,ii ). (3) 



1 da( L °) 



2D q ^(z,»p) + ^pW(z)log 



<7 dz 

= 2D q ^(z,iJ,p) + C^(z,fXp). 
The . . . represent terms which are well behaved as z — > 1. 
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In this equation, all unknown long-distance effects are related to the behaviour of D q ^(z, /j,q), 
the initial value of this fragmentation function which has been fitted to the data at some 
initial scale /i in [[|. As D q _,^(z, //p) is exact, this solution does not take the commonly 
implemented [§, || resummations of \og(fi F ) into account and when used to evaluate the 
photon + 1 jet rate at 0(a) yields a factorization scale independent prediction for the cross 
section. 

In the Durham jet algorithm, s cut ~ sz(l — z) 2 / (1 + z) ~ p\ |7j where pt is the transverse 
momentum of the photon with respect to the cluster. For z < 1, we find that fi 2 F ~ 
s cut and /if* ^> //q ■ The 'direct' contribution in eq. (g) is therefore suppressed relative 
to the fragmentation contribution. The conventional assignment of a power of l/a s to 
the fragmentation function can in this case be motivated, this contribution is indeed more 
significant. However, as z —>■ 1, we see that the transverse size of the photon jet cluster 
decreases such that s cut — > 0. The hierarchy s cut ~ fMp and \j? F 3> //q is no longer preserved 
and both contributions in eq. (Q) are important. Large logarithms of (1 — z) become the 
most dominant contributions. Being primarily interested in the high z region, in [Q] it was 
chosen not to impose the conventional prejudice and resum the logarithms of fi F a priori 
but to work within a fixed order framework, to isolate the relevant large logarithms. 

We have performed the calculation of the 0(a s ) corrections to the 'photon' + 1 jet rate 
using the same democratic procedure to define the photon as in The details of this 

fixed order calculation have been presented in ||. In the following, we shall limit ourselves 
to outline the main characteristics of this calculation, to summarize the results and to show 
how these compare with the available experimental data from ALEPH. 

2.2 The 'photon' + 1 jet rate at 0(aa s ) 

The 'photon' + 1 jet rate in e + e~ annihilation at 0(aa s ) receives contributions from five 
parton-level subprocesses displayed in Fig. [[]. Although the 'photon' + 1 jet cross section is 
finite at 0(aa s ), all these contributions contain divergences (when the photon and/or the 
gluon are collinear with one of the quarks, when the gluon is soft or since the bare quark-to- 
photon fragmentation function contains infinite counter terms). All these divergences have 
to be isolated and cancelled analytically before the 'photon' + 1 jet cross section can be 
evaluated numerically. 

The various configurations where the tree level process 7* — > qqg^ contributes to the pho- 
ton + 1 jet rate are illustrated in Fig. The associated real contributions can be separated 
into three categories: either theoretically resolved, single unresolved or double unresolved. 
Within each singular region which is defined using a theoretical criterion s m ; n , the matrix 
elements are approximated and the unresolved variables analytically integrated. The eval- 
uation of the singular contributions associated with the process 7* — > qqg'j is of particular 
interest as it contains various ingredients which could directly be applied to the calculation 
of jet observables at next-to-next-to-leading order. Indeed, besides the contributions arising 
when one final state gluon is collinear or soft (single unresolved contributions, see fig.||(b)), 
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Figure 1: Parton level subprocesses contributing to the photon + 1 jet rate at 0(aa s ). 



there are also contributions where two of the final state partons are theoretically unresolved, 
see fig.|2|(c) . The three different double unresolved contributions which occur in this cal- 
culation are: the triple collinear contributions, arising when the photon and the gluon are 
simultaneously collinear to one of the quarks, the soft/ collinear contributions arising when 
the photon is collinear to one of the quarks while the gluon is soft and the double single 
collinear contributions, resulting when the photon is collinear to one of the quarks while the 
gluon is collinear to the other. A detailed derivation of each of these singular real contribu- 
tions and of the singular contributions arising in the processes depicted in Fig. |TJ(b)-(d) has 
been presented in [||. 

Combining all unresolved contributions present in the processes shown in Fig. |l|(a)-(d) 
yields a result that still contains single and double poles in e. These pole terms are however 
proportional to the universal next-to-leading order splitting function P^> || and a convolu- 
tion of two lowest order splitting functions, {Pgq ®P$)- Hence, they can be factorized into 
the next-to-leading order counterterm of the bare quark-to photon fragmentation function 
PTJfl present in the contribution depicted in Fig. |I|(e), yielding a finite and factorization scale 
(Hf) dependent result ||. 

We have then chosen to evaluate the remaining finite contributions numerically using the 
hybrid subtraction method, a generalization of the phase space slicing procedure |Tl|, [12[]. The 
latter procedure turns out to be inappropriate when more than one particle is potentially 
unresolved. Indeed, in our calculation we found areas in the four parton phase space which 
belong simultaneously to two different single collinear regions. Those areas cannot be treated 
correctly within the phase space slicing procedure. Within the hybrid subtraction method 
developed in [|T3[| , a parton resolution criterion s min is used to separate the phase space 



into different resolved and unresolved regions. Phase space slicing and hybrid subtraction 
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Figure 2: Different final state 'photon' + 1 jet topologies arising from the tree level 7* — ► qq'yg 
process. The 'photon' jet is moving to the left while the recoiling hadronic jet moves to the 
right. Square brackets denote theoretically unresolved particles, round brackets experimental 
clusters. 



methods vary only in the numerical treatment of the unresolved regions. While the matrix 
elements are set to zero in the former method, one considers the difference between the 
full matrix element and its approximation in all unresolved regions in the latter. The non- 
singular contributions are calculated using Monte Carlo methods like in the phase space 
slicing scheme. 

The numerical program finally evaluating the 'photon' + 1 jet rate at 0(aa s ) contains 
four separate contributions. Each of them depends logarithmically (in fact as log 3 (y min )) on 
the theoretical resolution parameter y min = s m i n /Q 2 . The physical 'photon' + 1 jet cross 
section, which is the sum of all four contributions, must of course be independent of the 
choice of y m i n , the latter being just an artefact of the theoretical calculation. In Fig. || we 
see that the cross section approaches (within numerical errors) a constant value provided 
that 2/ m i n is chosen small enough, indicating that a complete cancellation of all powers of 
log(2/min) takes place. This provides a strong check on the correctness of our results and on 
the consistency of our approach. 

Finally, after factorization of the quark-photon singularities, the 0(aa s ) cross section 
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Figure 3: 0(aa s ) individual contributions (left) and sum of all 0(aa s ) contributions (right) 
to the photon + 1 jet rate for a single quark of charge e q such that ae q = 2tt, a s (N 2 — l)/2N = 
2n using the Durham jet algorithm with y cut = 0.1, and integrated for z > 0.7. 



takes the following form, 

1 dcr^ L °) _ 1 d(x( LO ) 
o"o dz 




CW(z,ti F ) + Cf <g> D q ^{z^ F ). (4) 



The lowest order cross section has been given in eq.(0) while the hard scattering coefficient 
functions C 4 - appearing explicitly in this equation are defined as follows. The (finite) next- 
to-leading order coefficient function is obtained numerically after the next-to-leading 
quark-photon singularity has been subtracted. More precisely, is obtained after summing 
all contributions which are independent of D q _> 7 (z, /j, f ) arising from the Feynman diagrams 
depicted in Fig. [I] together. A detailed description of the evaluation of has been given 
in H]. The coefficient function is the finite part associated with the sum of real and 
virtual gluon contributions to the process e + e~ — > qq. It is straightforward to evaluate, and 
can be found for example in [p 



2.3 Comparison with Experimental Data 

A comparison between the measured 'photon' + 1 jet rate and our calculation yielded 
a first determination of the quark-to-photon fragmentation function accurate up to 0(aa s ) 



15| . This function, which parameterizes the perturbatively incalculable long-distance effects, 
has to satisfy a perturbative evolution equation in the factorization scale fi F . At next-to- 
leading order (0(aa s )) this equation reads, 
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Pq^l q and Pq-ty are respectively the lowest order quark-to-quark and the next-to-leading 
order quark-to-photon universal splitting functions || [16], |17| . The next-to- leading order 
fragmentation function can be expressed as an exact solution of this evolution equation up 
to 0(aa s ) §, 

DM _ ^, wlog (j) + ^|(^_i)^ wlog (j) 

log g) C(l) ^^» wlog (|) 

2 \ 





a. 




V 2iV 


a. 


1 


2vr 


l 2N 



\og[^\ P^(z)®D(z,^) +D(z,fM ). (6) 



..2 " 99 



The initial function D^°\z,fi ) has been fitted to the ALEPH 1 jet data pf for ig, 



for the jet resolution parameter y cut = 0.06 and a s (M|) = 0.124 to yield Q, 

D^O^o) = |f (-P#>(z)log(l - zf + 20.8 (1 - z) - 11.07) , (7) 
where /zq = 0.64 GeV. The next-to-leading order (MS) quark-to-photon fragmentation func- 



tion (for a quark of unit charge) at a factorization scale ftp = Mz were shown in |15j and 
compared with the lowest order fragmentation function obtained in ||. A large difference 
between the leading and next-to-leading order quark-to-photon fragmentation functions was 
observed only for z close to 1, indicating the presence of large log(l — z) terms. 

Moreover, a comparison between the ALEPH data and the results of the 0(aa s ) calcula- 
tion using the fitted next-to-leading order fragmentation function for different values of y cu t 



can be found in || [H| . The next-to-leading order corrections were found to be moderate for 
all values of y cut demonstrating the perturbative stability of our fixed order approach. To 
test the generality of our results, we have considered two further applications: the 'isolated' 
photon rate and the inclusive photon distribution which we shall now briefly present. 

Using the results of the calculation of the photon + 1 jet rate at 0(aa s ) and the fitted 
quark-to-photon fragmentation function, we have determined the isolated rate defined as the 
photon + 1 jet rate for z > 0.95 in the democratic approach. The result of this calculation 
compared with data from ALEPH || and the leading order calculation || is shown in 
Fig. |j. It can clearly be seen that inclusion of the next-to-leading order corrections improves 
the agreement between data and theory over the whole range of y cut . It is also apparent 
that the next-to-leading order corrections to the isolated photon + 1 jet rate obtained in 
this democratic clustering approach are of reasonable size indicating a good perturbative 
stability of this isolated photon definition. 

The OPAL collaboration has recently measured the inclusive photon distribution for final 
state photons with energies between 10 and 42 GeV 0]. Fig. ^| shows our (scale independent) 
predictions for the inclusive photon energy distribution at both leading and next-to-leading 

1 Note that the logarithmic term proportional to P$ (z) is introduced to ensure that the predicted z 
distribution is well behaved as z — -> 1 101. 
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Figure 4: The integrated photon + 1 jet rate above z = 0.95 as function of y cn t, compared 
with the 0(a) and 0(aa s ) order calculations including the appropriate quark-to-photon 
fragmentation functions. 



order. We see good agreement with the data, even though the phase space relevant for 
the OPAL data, which corresponds to x 7 values as small as 0.2, far exceeds that used to 
determine the fragmentation functions from the ALEPH photon + 1 jet data. 



3 The photon + 1 jet rate in the conventional approach 

In the conventional approach, the parton-to-photon fragmentation function -Dj_> 7 satisfies an 
all order inhomogeneous evolution equation JT7| . Usually, these equations can be diagonalized 
in terms of the singlet and non-singlet quark fragmentation functions as well as the gluon 
fragmentation function. However when analyzing the global features of the solutions of these 
evolution equations, as was discussed in [pi|, several simplifications can be consistently made. 
For example, the gluon-to-photon fragmentation function is by orders of magnitude smaller 



than the quark-to-photon fragmentation functions, as was shown in |18[ . Its contribution to 
the photon cross section can safely be ignored. Consequently, the flavour singlet and non- 
singlet quark-to-photon fragmentation functions become equal to a unique fragmentation 
function D q ^ which satisfies an evolution equation having a similar form than the next- 
to- leading order evolution valid in the fixed order approach, see eq. Unlike in eq. (|5|) 
though, the strong coupling a s is now a function of the factorization scale, it runs. 
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Figure 5: The inclusive photon energy distribution normalized to the hadronic cross section 
as measured by the OPAL Collaboration compared with the 0(a) and (D(aa s ) order calcu- 
lations including the appropriate quark-to-photon fragmentation functions determined using 
the ALEPH photon + 1 jet data. 



The full solution D q ^ of the inhomogeneous evolution equation is given by the sum 
of two contributions; a pointlike (or perturbative) part -Dqi» 7 which is a solution of the 
inhomogeneous equation follows eq. (H) and a hadronic (or non-perturbative) part which 
is the solution of the corresponding homogeneous equation. In the conventional approach, 
approximate solutions of these evolution equations are commonly obtained as follows || || . 
First an analytic solution in moment space is obtained in the leading logarithm (LL) or 
beyond leading logarithm (BLL) approximations. These are then inverted numerically to 
give the fragmentation function in x-space. At LL only terms of the form (a™ ln™ +1 jjb 2 F ) are 
kept while at BLL both leading (a™ ln n+ jj? F ) and subleading (a™ In™ fi F ) logarithms of the 
mass factorization scale \ip are resummed to all orders in the strong coupling a s . 

It is worth noting that both LL and BLL solutions have an asymptotic behaviour given 

by ' 

where a(z) contains the splitting function -Pj_L- This asymptotic form lends support to the 
common assumption that the quark-to-photon fragmentation function D q ^ is 0(a/a s ). 
This assumption is in contrast with that adopted in the fixed order approach (cf. Section 0) 
where the quark-to-photon fragmentation function is 0(a). It leads to significant differences 
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in the respective expressions of the one-photon production cross sections. Indeed, the LL and 
BLL expression of the cross section in the MS scheme arising when one uses the corresponding 
resummed LL or BLL fragmentation functions in this approach reads 

1 da LL 



<j dz 

= + (g) C<p<8>Dm(z,p F ) + MJ CV>(z,i* F ). (9) 

As can be seen from these equations, no direct term contributes to the cross section at the 
LL level, while only the 0{a) direct term contributes at the BLL level to it. At the 
BLL level, contributions arising from|l|(a)-(b) do not enter in the cross section. As explained 
at length in [ 15(1 ) this conventional procedure of associating an inverse power of a s with the 



fragmentation function is clearly appropriate when the logarithms of the factorization scale 
fiF are the only potentially large logarithms but is problematic when different classes of large 
logarithms can occur as it is the case in the 'photon' + 1 jet cross section. 

All unknown long-distance effects are related to the behaviour of the input fragmentation 
function DgP(z,fJ>o) implicitly present in eq.(^). In the approaches of GRV or BFG, the 
non perturbative input function D™p(z, /xq) is treated with only minor differences. Those 



have been detailed in fllgfl . We shall here concentrate on the major common points in 
these approaches. At LL both GRV and BFG agree that D™p(z, /Xq) is negligible and 
can be described by a vector meson dominance model (VMD) as explained in 0] and |J 
respectively. However at BLL and in the MS scheme, the input fragmentation function 
cannot be negligible due to the presence of the direct term and cannot be described by 
a VMD input alone. Indeed, C^\z) diverges as z — > 1 and would drive the cross section to 
unacceptable negative values if a VMD input alone is considered for the input fragmentation 
function. Note that the requirement that the cross section is positive led the authors in 
0, T3| to consider in the fixed order approach a term proportional to P^ ln(l — z) 2 in the 
expression of D^^z, /Jo)- In summary, in any resummed or fixed order approach, as soon as 
the direct term enters the cross section, as it does in the MS factorization scheme, the 
input fragmentation function Dq?(z, /io) must compensate the large z behaviour of C^°\ 
Consequently, this input fragmentation function D^P(z,fio) as well as the total solution 
D q ^ 1 (z, /ip) in this MS scheme should clearly exhibit a divergent behaviour as z — > 1 . 

In Fig. |5] we compare the analytic expression of the fragmentation function obtained in 
the fixed order approach, eq.(^) with the BLL parameterizations of GRV and BFG for the 
numerically resummed solutions. We clearly see, that the fixed order solution does diverge 
as z — > 1 while the numerical solutions do not. This significant disagreement is mainly due 
to deficiencies in the numerical parameterizations. In fact, it can be traced back to the 
presence of logarithms of (1 — z) that are explicit in the expanded result. These logarithms 
should also be present in the numerical resummed results. However, the parameterizations 
are necessarily obtained by inverting only a finite number of moments and it is a well known 
problem to describe a logarithmic behaviour with a polynomial expansion. This clearly 
indicates that the presently available parameterizations for the resummed fragmentation 
functions are not accurate at large z and particularly for z > 0.95. 
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Figure 6: The quark-to-photon fragmentation function zD u _, 7 (z, fip) evaluated at /ip = M z 
in the (MS)-scheme. The NLO fit from the ALEPH 'photon' + 1 jet data is shown as solid 
line. The GRV (BFG) parameterization is shown dashed (dotted). 



Except in the very high z region however, we see that, the various fragmentation functions 
generally agree well with each other in shape and magnitude. Consequently we can expect, 
that predictions for the inclusive photon cross sections (which run over a wide range of z) 
will be largely in agreement, while significant differences may be apparent in the 'photon' + 
1 jet estimates which focus on the large z region. Indeed, we mentioned that the OPAL data 
were in agreement with predictions using either (GRV or BFG) parameterizations in the 
conventional approach, and we showed in Section [| that these data were also in agreement 
with the predictions obtained in our fixed order approach. 

Let us now concentrate on the 'photon' + 1 jet cross section, an observable which is 
sensitive on the large z region (0.7 < z < 1) and see how the predictions obtained in a 
conventional approach compare with the ALEPH data. In the following, we focus on one 
particular value of the jet clustering parameter y cut , y cut = 0.1. Fig. |7| shows the BLL 
'photon' + 1 jet rates obtained using either the GRV or BFG parameterizations of the 
photon fragmentation function for fip = Mz- Ignoring the large z region where we have 
reason to doubt the accuracy, we see that the BFG predictions lie systematically below that 
obtained using the GRV parameterization and go through the experimental data points. As 
discussed in |18| , this difference is due to both the choice of hadronic scale and the non-VMD 



input. The BFG input is smaller and the 'photon' + 1 jet data clearly selects this choice. 
Notice however, that the BFG parameterization for the fragmentation function unlike that 
of the GRV group was proposed well after the ALEPH data were released. 
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Figure 7: The 'photon' + 1 jet rate in the conventional approach using the GRV or the BFG 
inputs. The experimental data is taken from ||. 



4 Conclusions 



In summary, in Section ^| we have outlined the main features of the calculation || of the 
'photon' + 1 jet rate at 0(aa s ). Although only next-to-leading order in perturbation theory, 
this calculation contains several ingredients appropriate to the calculation of jet observables 
at next-to-next-to-leading order. In particular, it requires to generalize the phase space slic- 



ing method of |TT] , |l2f to take into account contributions where more than one theoretically 
unresolved particle is present in the final state. The 'photon' + 1 jet rate has then been eval- 
uated for a democratic clustering algorithm with a Monte Carlo program using the hybrid 



subtraction method of |L3| . The results of our calculation, when compared to the data || on 
the 'photon' + 1 jet rate obtained by ALEPH, enabled a first determination of the process 
independent quark-to-photon fragmentation function at 0(aa s ) in a fixed order approach. 
As a first application, we have used this function to calculate the 'isolated' photon + 1 jet 
rate in a democratic clustering approach at next-to-leading order. The inclusion of the QCD 
corrections improves the agreement between theoretical prediction and experimental data. 
Moreover, it was shown that these corrections are moderate, demonstrating the perturbative 
stability of this particular isolated photon definition. We have also shown that the inclusive 
photon energy distribution computed in this fixed order approach and using the quark-to- 
photon fragmentation function determined with the ALEPH data and is in good agreement 
with the recent OPAL data. 

In Section |3| we have outlined the main characteristics of the conventional approach have 
described how the LL and BLL solutions of the evolution equation are obtained. An impor- 



12 



tant feature characterizing this approach concerns the power of a s associated to the photon 
fragmentation function. Although this appears to be 0(a), inspection of the evolution equa- 
tion suggests a logarithmic growth with /i F , and in conventional approaches, it is ascribed a 
nominal power of a/a s . Consequences for the LL and BLL cross sections were also discussed 
in this section. 

The full /i^-dependent solution D q ^ n (z, fXp) of the evolution can only be obtained pro- 
vided some non-perturbative input is given. In the approaches of GRV or BFG, this input 
has two pieces, a small vector meson dominance contribution together with a perturbative 
counterterm. In either case, we have found that the large z behaviour of the fragmentation 
functions is not well reproduced by the parameterizations, the main problem being to de- 
scribe a logarithmic behaviour with a polynomial. Predictions using any of the presently 
available (GRV or BFG) resummed fragmentation functions do therefore not yield accurate 
results in the region z > 0.95. Ignoring this high z region, the BLL predictions for the 
photon + 1 jet rate obtained using the BFG parameterization are found to be in agreement 
with the ALEPH data, while the predictions using the fragmentation function of GRV lie 
too high. 

To summarize, we have seen that the inclusive and 'photon' + 1 jet data from LEP can 
be described using either the 0(aa s ) fragmentation function whose non-perturbative input 
is fitted to the ALEPH data or using the BLL parameterization of BFG. In the latter case, 
the agreement needs however to be restricted to z- values below 0.95. 
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